Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/5050--03-8655 


Antenna  Coupling/Isolation  Analysis: 
Generation  to  3-D 

L.  N.  Medgyesi-Mitschang 

Signature  Technology  Office 
Systems  Directorate 


February  20, 2003 

20030320  050 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 

Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

February  20, 2002  NRL  Memorandum  Report 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Antenna  Coupling/Isolation  Analysis:  Generation  to  3-D 

5b.  GRANT  NUMBER 

64-K267-L-3 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

L.N.  Medgyesi-Mitschang 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

Naval  Research  Laboratory 

4555  Overlook  Avenue,  SW 

Washington,  DC  20375-5320 

NRL/MR/5050-- 03-8655 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  VA  22217-5660 

11.  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  report  summarizes  a  generalized  formulation  for  three-dimensional  antenna  coupling  and  isolation  problems.  The  formulation  is  devel¬ 
oped  for  a  generic  transmit/receive  configuration.  The  formulation  is  based  on  representing  the  boundary  conditions  arising  in  the  problem  by 
using  integro-differential  representations  of  Maxwell’s  equations.  To  streamline  the  derivations,  integral  operators  are  introduced.  The  resulting 
system  of  coupled  equations  for  the  electric  and  magnetic  fields  is  solved  using  the  Method  of  Moments  (MM)  technique.  Specifically,  the 
Galerkin  variant  of  this  technique  is  adopted  here,  resulting  in  a  symmetric  system  matrix  for  the  problem. 

The  formulation  discussed  here  is  specialized  to  the  two-domain  case.  This  refers  to  problems  where  the  antenna  or  array  is  embedded  in  a 
conducting  surface  constituting  one  domain.  The  second  domain  is  an  intervening  region  that  is  penetrable.  Depending  on  the  nature  of  the 
penetrable  region,  it  provides  the  isolation  or  coupling  between  the  transmit  or  receive  antennas  and  arrays. 

The  resulting  generalized  formulation  has  been  implemented  with  computer  software  that  is  an  adaptation  of  the  MM-based  CARLOS  code. 
Discussion  of  the  associated  numerical  implementation  and  software  are  documented  elsewhere  and  omitted  from  this  report. 

15.  SUBJECT  TERMS 

Theoretical  formulation  for  3-D  antenna  analysis;  Method  of  Moments;  Antenna  coupling 

1  16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  ABSTRACT 

OF  PAGES 

L.N.  Medgyesi-Mitschang 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

UL 

16 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

(202)  404-3081 

Standard  Form  298  (Rev.  8-98) 

i 

Prescribed  by  ANSI  Std.  Z39.18 

CONTENTS 


1.  BACKGROUND .  1 

2.  STATEMENT  OF  PROBLEM . 2 

3.  3-D  FORMULATION  OF  THE  TWO-DOMAIN  PROBLEM .  5 

4.  SOLUTION  OF  THE  TWO-DOMAIN  PROBLEM . .  9 

5.  NUMERICAL  IMPLEMENTATION . 12 

REFERENCES . . . . . . .  13 

LIST  OF  FIGURES .  13 


ANTENNA  COUPLING/ISOLATION  ANALYSIS: 
GENERATION  TO  3-D 


1.  BACKGROUND 

In  an  earlier  report  [1],  we  discussed  the  issue  of  antenna  coupling/isolation  in  the 
presence  of  layered  materials  using  a  two  dimensional  (2-D)  theoretical  formulation. 
Such  a  formulation  can  provide  general  trends  for  antenna  isolation.  In  particular,  one  can 
study  effects  of  discontinuities  and  tapers  of  the  material  layers  used  for  isolation. 
Because  the  formulation  was  2-D,  the  numerical  implementation  results  in  algorithms 
that  are  inherently  computationally  fast  and  require  minimum  computer  memory. 

After  general  trends  are  established  with  a  2-D  simulation,  a  more  accurate  calculation 
requires  a  3-D  analysis.  Such  an  analysis  is  important  when  cross  polarization  effects  are 
present.  Unfortunately,  3-D  simulations  are  much  slower  and  require  substantially  more 
computer  resources.  The  following  summarizes  the  generalization  of  the  2-D  formulation 
to  the  3-D  case.  For  clarity  and  ease  of  reference,  we  adhere  to  the  notational  scheme  in 
[1]  as  much  as  possible. 
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2.  STATEMENT  OF  PROBLEM 


Let  us  consider  a  radiating  and  a  receiving  system  embedded  in  a  surface  as  depicted  in 
Figure  1.  Again  for  generality  each  of  these  systems  may  consist  of  individual  or  multiple 
antenna  elements.  The  surface  between  the  elements  can  be  conducting  or  consist  of 
single-  or  multi-layered  dielectric  (and/or  magnetic)  coatings  on  a  conducting  surface. 
The  electrical  characteristics  of  this  surface  strongly  dictate  the  isolation  (or  coupling)  of 
the  transmitting  and  receiving  antenna  elements. 

As  in  the  2-D  formulation,  the  outer  surface  of  these  coatings  can  be  smooth  or 
discontinuous  including  corrugated  conducting  surfaces.  Such  a  surface  provides  a 
reactive  boundary  that  leads  to  surface  wave  guidance.  With  an  appropriate  choice  of 
profiles,  surfaces  can  be  designed  to  have  stop-band  or  pass-band  properties.  The 
antenna/isolation  material  regions  can  be  flat  or  curved.  In  most  cases  in  practice  the 
array  surface  is  flat.  The  intervening  region  between  the  transmitting  or  receiving  array 
may  be  flat  or  curved.  ° 

To  provide  a  framework  for  the  subsequent  discussion,  we  depict  the  antenna  isolation 
problem  in  Figure  2.  The  discontinuities  and  terminations  associated  with  the  antenna 
apertures  and  transitions  between  the  coatings  can  be  either  electrically  large  or  small 
with  respect  to  wavelength.  The  transmitting  and  receiving  antennas/arrays  may  have 
different  footprints  as  well. 

The  problem  depicted  in  Figure  2  is  sometimes  called  a  two-domain  problem.  The 
antenna  or  array,  embedded  in  the  conducting  surface,  is  designated  as  one  domain  (or 
region).  Here  it  is  with  bounded  by  the  surfaces  S0  and  S2.  The  dielectric  region  is 
designated  as  the  second  domain  (or  region)  bounded  by  the  surface  S2and  with  the 
surface  S,  at  the  free  space/dielectric  interface.  Clearly,  more  complex  arrangements  are 
possible.  These  can  also  be  treated  with  the  following  analysis.  They  do  not  introduce  any 
new  theoretical  or  mathematical  issues  other  than  more  notational  complexity  in  the 
derivations.  In  the  interest  of  clarity  and  simplicity  we  will  forego  such  generalizations. 
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Conducting  enclosure 


Region  Rq 

Figure  1 :  Generic  Transmit/Receive  Configuration 
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R\  (free  space  region  : 
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Figure  2:  Detail  of  Two-Region  Problem 
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3. 3-D  FORMULATION  OF  THE  TWO-DOMAIN  PROBLEM 


Referring  to  Figure  2,  the  transmitting  and  receiving  antennas  or  arrays  are  embedded  in 
the  conducting  surface  S„.  The  conducting  region  R0  need  not  be  a  solid  conductor.  Often 
in  practice  this  is  an  enclosure  containing  the  antenna  connectors,  phase  shifters,  and 
other  electronic  components  of  the  array  back-plane.  As  in  the  earlier  2-D  formulation, 
we  assume  that  all  the  antenna/array  elements  themselves  are  surface  embedded.  The  LID 
program  assumes  such  elements. 

Region  R2  defines  the  penetrable  region  with  the  exterior  free  space  boundary  of  Sj  and 
the  remaining  boundary  being  surface  S2.  The  dielectric  region  is  characterized  by 
complex  permittivity  and/or  permeability.  For  the  moment  to  avoid  unnecessary 
notational  complexity,  we  assume  that  this  region  is  a  homogeneous,  isotropic  layer.  The 
formulation  can  be,  and  has  been,  generalized  for  the  case  of  multi-layered  and  tapered 
media.  The  CARLOS  computer  implementation  incorporates  these  generalizations. 

Using  the  notation  in  Figure  2,  the  defining  equations  for  the  electromagnetic  interactions 
for  each  of  the  regions  can  be  written  in  terms  of  integral  equations  written  in  integral 
operator  form.  Topologically  they  have  the  same  form  as  the  2-D  formulation  discussed 
earlier  in  [1].  But  in  this  case  the  spatial  coordinate  becomes  r.  The  total  electric  and 
magnetic  fields  in  the  free-space  region  R,  are  given  respectively  by: 

ex(r)E,(r)  =  -L\(JX  +  Jxd+)(r)  +  KXMX  (r),  reRx  (1) 

dx(r)Hx(r)  =  -KX(JX  +  Jxd+)(r) -Ll.MUU  r  €  Rx  (2) 

Vi 

The  vector  quantities  or  functionals  J(r)  and  M(r)  refer  to  the  electric  and  magnetic 
surface  currents  on  the  respective  boundaries  that  separate  the  respective  regions, 
depicted  in  Figure  2.  Fundamentally  these  are  the  Kirchhoff  loop  currents  of  circuit 
theory,  generalized  here  to  the  vector  field  problem.  The  currents  have  a  spatial 
dependence  r  where  r  is  the  spherical  coordinate  for  a  3-D  geometry. 

The  subscripts  on  the  current  terms  refer  to  the  surfaces  forming  the  various  boundaries. 
The  superscripts  refer  to  the  type  of  the  boundary  (i.e.,  conducting  or  dielectric)  on  which 
the  currents  reside.  For  example,  Jx  refers  to  the  electric  currents  on  a  conducting 
surface.  Likewise,  Jf+  refers  to  the  currents  on  the  upper  surface  of  the  dielectric  region 
forming  the  boundary  with  the  free  space  region. 

In  like  manner  the  total  fields  in  the  penetrable  region  R2  can  be  written  as: 

62(r)E2(r)  =  -El(~J2  +  Jdx~)(r)  +  K2Mx(7),  r  e  R2  (3) 
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The  symbol  0(r)  denotes  the  usual  Heaviside  function  on  and  inside  the  boundary  of 
region  R,  or  R2.  For  region  R,  it  is  defined  as: 


0x(r)  = 


1 

•  1/2 
0 


re/?, 
redRx  • 
otherwise 


(5) 


To  achieve  notational  compactness  we  have  introduced  the  integro-differential  operators 
L  and  K  in  Eqs.  (1-4),  defined  as 


LtXfr)  =  J  (j(O^X(r)  +  -j-  VV'  •  X(r))^  (r  -  r')dr' 
as,  wei 


(6) 


KxX(r)  =  J  X(r')  x  V<J>,(r  -  f')df' 
as, 


(7) 


where  the  vector  functional  X(r)  refers  to  either  the  electric  J(r)  or  the  magnetic  M(r) 
currents,  respectively. 


For  the  3-D  case,  the  Green  s  function  for  region  Rj  (i.e.,  the  free-space  region)  is: 

=  7'\  l exP[~iW ~ (8) 

where  k,  is  the  wave  number  in  region  R,. 

For  the  penetrable  (dielectric)  region  R2 ,  the  corresponding  Green's  function  is  defined  in 
terms  of  the  permittivity  and/or  permeability  and  the  wave  number  of  that  region.  The 
integrals  are  over  the  surfaces  or  boundaries  between  the  regions.  The  integrals  in  the  L 
and  K  operators  are  defined  in  the  Cauchy  principal  value  sense. 

The  boundary  conditions  on  each  side  of  the  penetrable  surface  requires  the  continuity  of 
the  tangential  fields,  namely: 


E^hanS, 

=  E2  (r)tanS, 

(9) 

WLs, 

=  ”2  (r)lanS> 

(10) 
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The  tangential  electric  fields  vanish  on  the  conducting  surfaces  yielding 


WLs,=°  do 

£,<r)ms,=  0  (12) 


Letting  =  Mj+  =  -Aff  and  Jf  =  Jf+  =-J\  on  the  dielectric/free-space  boundary, 
one  obtains  six  coupled  integral  equations  for  the  four  unknown  currents  on  all  the 
boundaries.  For  points  on  the  dielectric/free-space  boundary,  we  obtain 

\hJcx(r)  +  (L,+  L, )j}(. r )  -  ^/2( r )  -  (Kx  +  K^MJr)]^  =  0  (13) 


KlJ1c(r)  +  (K1+K2)jf(r)-K2J2(r)-(\Ll  +\L2)Mx(r) 

Vi  Vi 


=  0 


(14) 


I  tan  Sl 


On  the  exterior  conducting  boundary,  electric  and  magnetic  field  integral  equations  (EFIE 
and  MFIE)  are  respectively, 


[Ll(Jlc  +  J1d)(r)-R 'A(r)]tanS  =0 


(15) 


-  +  «j  x  KX(JX  +  J? )(r)  +  -4 nx  x  LyMfr) 
2  Vi 


=  0 


(16) 


tanSn 


Finally,  on  the  interior  conducting  boundary,  the  EFIE  and  MFIE  formulations  are 
respectively. 


[L2(j2-J?)(r)  +  K2Mx(r)\anS  =0 


(17) 


^  *  KiCh  ~  Jf)(r)  - X  h.Mx(r) 

2  v2 


=  0 


(18) 


l  tan  S-, 


The  imposition  of  the  boundary  conditions  results  in  six  equations  for  the  four  unknown 
currents  on  the  boundaries.  Subsequently  only  four  equations  will  be  used,  i.e.  those 
associated  with  the  PMCHW  formulation  (after  Poggio,  Miller,  Chu,  Harrington  and 
Wu),  namely  Eqs.  (13)-(14),  and  the  E-field  boundary  condition  on  the  conducting 
surfaces,  i.e.,  Eqs.  (15)  and  (17).  [2,3] 


7 
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configuration  was  primarily  chosen  to  introduce  the  notation  required  to  carry  out  the 
above  derivation. 

The  principal  restriction  so  far  is  that  the  problem  consists  of  two  domains:  one 
conducting  and  the  other  penetrable  and  homogeneous.  As  noted  before,  this  restriction 
was  imposed  to  streamline  the  foregoing  derivation.  A  more  general  derivation  is 
possible  that  not  removes  this  restriction.  The  CARLOS  computer  implementation 
involves  the  more  general  case,  i.e.,  where  Region  R2  is  multi-layered.  A  further 
generalization  allows  the  antenna/array  elements  to  be  embedded  within  the  penetrable 
(dielectric)  medium  itself.  To  achieve  RCS  suppression  of  array  apertures  often  requires 

specialized  antenna  windows.  In  such  cases  the  foregoing  generalized  analysis  option  is 
required. 
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4.  SOLUTION  OF  THE  TWO-DOMAIN  PROBLEM 


There  are  a  number  of  ways  of  solving  the  derived  system  of  equations.  In  this  report  we 
confine  ourselves  to  the  method  of  moments  (MM)  technique.  We  use  the  Galerkin 
variant  of  this  method.  This  approach  has  proven  to  be  computationally  robust  and  highly 
flexible  in  handling  a  great  variety  of  problems.  [4] 

Following  the  Galerkin  MM  technique,  one  discretizes  the  surfaces  in  Figure  2  with  flat 
triangular  facets.  The  currents  on  each  facet  are  expanded  in  terms  of  roof-top  function  s 
at  each  interior  edge  on  a  surface.  Half  roof-tops  are  used  for  junctions  between  two 
surfaces  intersecting  where  for  instance  a  conducting  and  a  dielectric  surface  meet.  These 
"junction"  functions  allow  the  Kirchhoff  continuity  condition  to  be  met  in  the  numerical 
implementation.  A  detailed  description  of  these  functions  and  the  associated  junction 
treatments  is  given  in  Section  2.6  of  [5].  Next  one  expands  the  unknown  currents  on  the 
conducting  and  dielectric  boundaries  via  the  set  of  expansion  functions.  Here  we  choose 
triangle  functions  of  the  Rao-Wilton-Glisson  (RWG)  type. 

The  expansion  of  the  electric  and  magnetic  surface  currents  on  each  facet  can  be  written 
explicitly.  For  example. 


as) 

n 

M(t)  =  -rloJ£Knfn(r)  (20) 

n 

where  fn(r)  are  roof-top  functions,  and  ln  and  Kn  are  the  unknown  coefficients  of  the 
expansions.  Since  the  foregoing  expansions  are  over  a  finite  number  of  terms,  the  current 
representations  are  approximate.  However,  in  practice  if  the  facet  surface  discretization  is 
on  the  order  of  A/10  per  linear  dimension,  convergent  results  are  obtained. 

The  current  expansions  are  substituted  into  the  integral  equations  and  the  Galerkin  testing 
procedure  is  carried  out  in  Eqs.  (13)-(15)  and  (17).  The  details  are  discussed  in  [5].  The 
Galerkin  procedure  yields  a  system  of  equations  linear  in  the  unknown  coefficients  of  the 
electric  and  magnetic  currents  In  and  Kn  on  the  surfaces.  This  system  is  customarily 
written  in  matrix  form  as: 


ZI  =  V  (21) 

The  Z  matrix  contains  all  the  electromagnetic  interactions  of  the  transmitting  and 
receiving  antennas  in  the  presence  of  the  penetrable  region  R2  between  them.  It  should  be 
noted  that  these  interactions  include  both  the  radiated  field  coupling  between  the  antennas 
as  well  as  the  surface  current  (or  conduction)  coupling.  For  closely  collocated 
transmitting  and  receiving  systems,  the  latter  is  dominant  and  depends  on  the  electrical 
characteristics  of  the  surface  (ground  plane)  between  them. 
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The  general  form  of  the  Z  matrix  is: 


f  LfZ'S,;*!)) 


Z  = 


{ +L(S1,S1;R2)J 
rK(Sx,Sx;Rx)+> 
^  K(S\  !  ^2)  ) 


US*,SX ;  R^) 

>  R\) 


'-K(SX,SX;RX) 
K~E(SX,SX ; Rl)J 
'%L(SX,SX;RX)  + 

-K(S0,SX;RX) 

KfS^S,^) 


L(St,Scx;Rx) 

K(Sf,%;Rx) 

L(S0,S0;  Rx) 
0 


L(SX  ,S2;  R} ) 

~K(SX  ,S2; R2) 
0 

l/S2,S2i  R2) 


(22) 


The  submatrix  blocks  of  Z  were  obtained  by  forming  the  vector  interproduct  of  the 
integral  operators  with  the  testing  function.  In  the  case  of  the  Galerkin  implementation  of 
the  MM  technique,  the  latter  is  the  complex  conjugate  of  the  expansion  functions  used  for 
the  electric  and  magnetic  currents.  Thus  for  example,  the  block  L(SX  S2;R)  is  obtained  by 

forming  the  interproduct  of  the  integral  operator  L  on  S2and  the  testing  function  on  S,. 
Further  details  of  this  procedure  are  given  in  [5]. 


The  column  vector  I  in  Eq.  (21)  represents  the  coefficients  of  unknown  currents  on  the 
conducting  and  dielectric  boundaries.  Specifically,  it  is  given  as: 


4 


1  = 


Afj 

Jl 


h 


(23) 


Note  subparts  of  this  column  vector  are  the  electric  and  magnetic  unknown  currents  on 
the  various  surfaces.  The  superscripts  designate  the  latter.  The  column  vector  V  in  Eq. 
(21)  denotes  the  aperture  excitation  of  the  transmitting  antenna,  specifically  the  electric 
and  magnetic  fields  in  the  transmitting  aperture  (or  apertures).  It  is  given  as: 


V  =  -\ 


E(SX) 

ms,) 

E(SX) 


0 


(24) 


The  elements  of  V  are  given  known  quantities.  For  example,  if  the  antenna  aperture  is 
embedded  in  the  conducting  surface  S0,  the  quantity  E(S0)  represents  the  electric  field 
excitation  (or  voltage  impressed)  on  the  aperture.  In  this  case  the  other  entries  of  the  array 
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are  zero.  Knowing  the  aperture  excitation,  the  foregoing  matrix  system  in  Eq.  (21)  can  be 
solved  for  the  unknown  currents  either  by  LU  decomposition  or  an  iterative  technique. 
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5.  NUMERICAL  IMPLEMENTATION 

In  add  non,  (he  software  allows  the  distribution  of  near  and  ultra^lr  fields  to 
rmppecUrUo  tblT"' pr0ximi,y  of  anlennas- If  dcsired,  the  latter  fields  can 

fiatrrfe 
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Figure  1:  Generic  Transmit/Receive  Configuration  (Two-Domain  Problem) 

Figure  2:  Detail  of  Two-Domain  Problem 
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